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Based on a systematic investigation of J/il;{tp') VP, where V and P stand for hght vector 
and pseudoscalar mesons, we identify the role played by the electromagnetic (EM) transitions and 
intermediate meson loop transitions, which are essential ingredients for understanding the J/ip and 
Tp' couplings to VP. We show that on the one hand, the EM transitions have relatively larger 
interferences in ip' — > pyr and K*K + c.c. as explicitly shown by vector meson dominance (VMD). 
On the other hand, the strong decay of ^' receives relatively larger destructive interferences from 
the intermediate meson loop transitions. By identifying these mechanisms in an overall study of 
J/^{ip') VP, we provide a coherent understanding of the so-called "pn puzzle". 

PACS numbers: 12.40.Vv, 13.20.Gd, 13.25.-k 



I. INTRODUCTION 

The decay channels J/ip{ip') VP, which are suppressed in QCD due to the violation of hadronic- 
helicity conservation have attracted much attention in the past few decades. According to this 
selection rule, one expects the ratio BR^ip' pTr)/BR{J/ip — > pi:) ~ (Mj/^/M^)^ ^ 0.35 [l|, which 
turns out to be much larger than the experimental data 2], BR{ijj' pn) / B R{J / ip pir) ~ (0.2±0.1)%. 
This significant discrepancy is known as the so-called "p7r puzzle" . An alternative expression of the "p7r 
puzzle" is via the ratios between J /if) and tp' annihilating into three gluons and a single direct photon: 

_ BRi:^};' ^hadrons) ^ BRj^' ^ e+e-) ^ 

~ BR{J/i}j hadrons) ~ BR{J/%P ^ e+e-) " ^ > 

which is empirically called "12% rule". The puzzle arises from the violation of the above empirical rule 
in exclusive channels such as pi: and K* K -\- c.c, where the branching ratio fractions are found to be 
orders-of-magnitude smaller than the approximate "12%". 

Since the first observation of such a large deviation by Mark-II Collaboration in 1983 [31 , many theo- 
retical explanations have been proposed to decipher this puzzle d, H, 0, H, H, [13, HHj fl^-lisl- fl3. fill [TgI. 
113, [3- They can be classified into three categories: i) J/ t/)- enhancement hypothesis, which attributes 
the small i?-value to the enhanced branching fraction of J /ij) decays; ii) ^/I'-suppression hypothesis, which 
attributes the small i?-value to the small branching ratio of ij)' decays; iii) and other hypotheses which 
do not simply belon g to the above two categories. Unfortunately, so far none of those solutions has been 
indisputably agreed [Tgl. [20j. 

In this proceeding, we report our recent efforts on understanding this issue. We shall identify i) the 
role play by EM transition in J/-0('i/'') VP in a VMD model, which has relatively large interferences 
in pTT and K* K -f c.c. channel; and ii) mechanisms which suppress the strong decay amplitudes for ipi' — > 
VP. We emphasize that our analysis is based on so far the-state-of-art experimental measurements of 
J/tljiip') ~^ 0- The systematic exposed in this study can provide some insights into the charmonium 
hadronic decays on a much more general ground. 
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FIG. 1: Schematic diagrams for J/ipixj}') j* ^ VP. 



II. STEP 1: EM TRANSITIONS IN THE VMD MODEL 



The importance of EM transitions in J/^p^ip') — > VP can be recognized by exphcit experimental 
observations. For instance, the branching ratios for isospin-violating decays, J/iAip') pr], pr]' , etc are 
compatible with those isospin-conserving channels such as uir], uirf , and (/iTy etc [2|. This is an evidence 
showing that the strong decay amplitudes become suppressed and have the same order of magnitude as 
the isospin-violating amplitudes, i.e. EM and strong isospin-violating transition. 

The role of the EM can be separately investigated due to the available experimental data for vector 
meson radiative decays, i.e. w, p, (j), K*, J/tP and i/j' 0. Moreover, precise measurements of vector 
meson decays into lepton pairs such as e"*"e~ are also available. This allows a well-constraint on the 
coupling constants required in the VMD model, and only leaves an overall form factor which takes care 
of the off-shell couplings, to be determined by experimental data. Three independent EM transitions for 
Vi^V2P are illustrated by Fig. [H 

Typical effective Lagrangian for the V^P coupling are: 

^^7P = ^e^.c^d'^V'^d^A^P (2) 

My 

where V'^{— p, uj, <j), J/ip, 4'' ■ ■ ■) and are the vector meson and EM field, respectively; Afy is the 
vector meson mass; e^uai3 is the anti-symmetric Levi-Civita tensor. 
The V^* coupling is described in VMD model, 

cy, = '-^y.A' , (3) 

where eMy/fy is a direct photon- vector- meson coupling in Feynman diagram language, and the isospin 
1 and component of the EM field are both included. 

The invariant transition amplitude for Vi — > 7* ^ V2P can thus be expressed as: 

Mem = Ma + Mb + Mc 

Jy2 My I ]yi My 2 JVlJV2 Mp J 

where gp-yj is the coupling for the neutral pseudoscalar meson decay to two photons; J- a, Th and Tc 
denote the form factor corrections to the transition of figure [TJ A monopole (MP) form factor is adopted 
here, 

with A = 0.542 ± 0.008 GeV and A = 0.577 ± 0.011 GeV determined by the isospin violated channels 
J/^{tp') prj, prf, wtt", and 07r° with a constructive (MP-C) or destructive phase (MP-D) between Fig.[T] 
(a) and (b), respectively. The form factors are introduced because we think that the non-perturbative 
QCD effects may play an important role in the transition at J/ip energy scale. 

The form factor J-c appearing in Eqs. Q can be determined in 7*7* scatterings. A commonly adopted 
form factor is 
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FIG. 2: Schematic diagrams for J/i/i ^ (^P via strong interaction, where the production of different components 
of the pseudoscalar P is demonstrated via (a): SOZI process; (b) DOZI process; and (c) glueball production. 
Similar processes apply to other VP channels as described in the text. 



where = My-^ and = are the squared four-momenta carried by the time-hke photons. We 

assume that the A is the same as in Eq. ([5]), thus, J-c = J'a^b- 

It should be noted that in Fig. [T] the direct application of V^P couplings extracted from experimental 
data will avoid uncertainties arising from a 7 — > y — s- VP treatment. Unknown energy-dependence of 
those couplings can then be absorbed into an overall form factor J^{q^) for which the cut-off energy is 
determined by fitting those isospin-violating decay branching ratios, i.e. J/^{ip') —^ prj, prj' , lutt and (jjir. 

In fact, one can learn more from the isospin-violating channels. If the EM transition is the dominant 
transition mechanism, one can expect that the 12% will be reasonably respected given that the J/ip and 
"0' wavefunctions are normal cc of {IS) and (25), respectively, and no significant interferences from other 
processes. As shown in Table HI one indeed sees that the 12% rule is satisfied though the experimental 
values still have large uncertainties. There might be contributions from strong isospin-violating transi- 
tions. However, the present experimental results suggest that their interferences with the EM transitions 
in the isospin-violating channel are relatively small. 



Decay channels R {%) Exp. data (%) 





8.97 


11.5 ±5.0 


pt] 


9.44 


23.5 ± 17.8 


LUTV 


9.01 


5.0 ± 1.8 




7.41 


< 62.5 



TABLE I: Branching ratio fractions of t/j' — > 7* ^ VP over J/t/j ^ 7* VP for those isospin-violating channels. 
Here, we only show results with MP-C form factor. The last column is extracted from the experimental date (31 ■ 



III. STEP 2: PARAMETERIZE THE STRONG DECAY TRANSITIONS 

For those isospin-conserved decays, i.e. J l'ip{ip') wq, cjr/, 077, 0?/, pi: and K* K -\- cc, the strong 
and EM decay processes are mixed. Recalling that the antisymmetric tensor form is the only coupling 
for VVP, we thus parameterize the strong decays in a way similar to Refs. [Ml, [H, [H, [13] . Some basic 
quantities can be defined via Fig. [2l the strength of non-strange singly OZI disconnected process 57/^; 
the parameter reflecting the SU(3) flavor symmetry breaking effects R, and the parameter r describing 
the relative strength between the DOZI and SOZI transitions. The expressions for the parameterized 
strong decay amplitudes are listed in Tabic HH 

In Ref. [25|, different treatments for glueball— mixing are investigated, which are denoted by Schemes 
I, II, III. Since the glueball components in rj and 77' are rather small, the mixing effects will not change 
out results on the strong decays of J/'ijj{ip') VP. Details about the parameter deflnitions and mixings 
can be found in Ref. [25'| . 

In order to take into account the size effects from the spatial wavefunctions of the initial and flnal-state 
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Decay channels 


Transition amplitude M = [Mi + M2 + Ms) 


H 






9j/i,(,l,')RW^rx2 + R{1 + r)y2 + Z2]J^{P) 






t 


9j/i,(r + 2r)x2 + V2Rry2 + V2z2]T{P) 



p -K 




p'^TT^ or p^TV'^ 


gj/i,{i,')J^{P) 


K*°K° or K*°K° 




K*+K- or K*-K+ 





TABLE II: General expressions for the transition amplitudes for J/ij}(^') — > VP via strong interactions. Param- 
eter Qj/^ and g^i axe proportional to the charmonium wavefunctions at origin and have different values for J/i/; 
and ^j)' , respectively. For 77 and rj' , a glueball mixing is also considered in the wavefunctions. 





ffj/^(xlO-^) 


ff^'(xl0-^) 


Scheme-I 


18.66 ±0.63 


2.22 ±0.49 


Scheme-II 


18.45 ± 0.70 


2.54 ±0.44 


Scheme-Ill 


17.52 ±0.66 


2.57 ±0.42 



TABLE III: Extracted coupling strengths for the SOZI transitions for three different 77 — 77'— glueball mixing 
schemes with a constructive mode for the EM amplitudes [25j . 

mesons, we apply the commonly used form factor 

:r2(P) = |P|2; cxp(-PV8/32), (7) 

where P and the I are the three momentum and the relative orbit angular momentum of the final- 
state mesons, respectively, in the J/^'i'ip') rest frame. We adopt f3 = 0.5 GeV, which is the same as 
Refs. [2i,[23,[2i,[2|. 

In Table IIIIl we list the values of the strong coupling strengths and g^' which are extracted 

by overall fittings to the isospin-conserved decay channels of J/'ip{ip') VP data including the EM 
transitions determined in the previous section. The predominant feature is that both values are stable 
and insensitive to the 77 — ry'— glueball mixing schemes. In case of the absence of the EM contributions, 
the "12% rule" fraction should be proportional to [gj/^/g^')'^- By taking the average of the squared 
values of Table Hill we obtain {gj/tfj/g^'Y — 1-8% which is much less than the expectation of the "12% 
rule", but larger than the experimental data, ^ (0.2 ± 0.1)%. 

This is not at all a trivial outcome. Several points can be learned here: i) The suppression to the ip' 
strong decay coupling is not exclusively on ij:' —> pn. Such a suppression is an overall effect on all the 
exclusive decays, ii) Due to the suppression on the strong decay coupling of the the EM transition 
amplitudes become compatible with the strong decay amplitudes with which the interferences produce 
deviations from naive expectations based on single transition mechanism. To be more specific, due to 
the interference, the pir decay is further suppressed, i.e. causes the so-called "/ott puzzle". The neutral 
K*'^K'^ _l_ (, jjg^g larger branching ratio than the charged one K*+K- + c.c. 0. iii) As shown in Fig. 
the DOZI transitions contribute to the isoscalar channels. This suggests that the exclusive decays have 
different features compared with the inclusive one from which the "12% rule" is embedded. 

In Table HVl we list the branching ratio fractions for those isospin-conserved channels and compare them 
with the experimental data. Within the experimental uncertainties, our results are in good agreement 
with the data. We also show the branching ratio fractions for exclusive EM transitions, and again, one 
can see that the "12% rule" is reasonably respected for exclusive transitions. 

In brief, the parametrization identifies the mechanism which not only causes puzzle in pi: channel, but 
also plays a role in other VP channels, since the EM transitions in pn and K* K + c.c. are relatively large 
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Decay channels R^''{%) Scheme I(%) Scheme II(%) Scheme III(%) Exp.data (%) 





8 


0.12 


0.16 


0.20 




pvr 


8 


0.12 


0.15 


0.19 


0.2 ±0.1 


ujr] 


8 


0.40 


0.35 


0.28 


< 0.6 ±0.1 


LOT] 


8 


5.33 


0.11 


0.29 


18.5 ± 13.2 


<i>v 


10 


2.78 


2.93 


3.30 


4.1 ± 1.6 




10 


5.00 


5.34 


8.86 


8.7 ±5.5 


K*+K- + c.c. 


9 


0.45 


0.41 


0.39 


0.4 ± 0.2 


K*°K° + c.c. 


9 


2.74 


2.79 


2.67 


2.7±0.7 



TABLE IV: Branching ratio fractions for ?/)' ^ 7* VP over J/tp — 7* VP for different rj-ri'— glueball with 
a MP-C form factor. R^^ denotes the ratios with exclusive EM transitions. The last column is extracted from 
the experimental date 01 . The stars in p^n^ channel denotes the unavailability of the data. 



due to large couplings for pnj and K* Kj, interferences between the suppressed strong decay am plit udes 
and enhanced EM amplitudes produce significant deviations from the expectation of "12% rule" [25l.[30l|. 

IV. STEP 3: MECHANISM SUPPRESSING THE STRONG DECAY AMPLITUDES OF 

■4,' VP 

Now, the last bit of the whole scenario comes to the point, "Why, and how the strong decay coupling 
g^' is suppressed?" In order to demonstrate this, we express the decay amplitudes as 



IvIjI,' 

which again benefits from the property of the antisymmetric tensor coupling among VVP fields. In the 
above equation, gj/^ and g^pt are real numbers fixed by Step 2 [25l|, while gjjlp and g'^T are the EM 
couplings fixed by Step 1 with relative phase angles and S^i fixed in Step 2. Detailed discussions on 
the phase angles can be found in Rcf. [25] , of which the values can be compared with those from Rcf. [3l| . 
Since any possible mechanism must contribute to the coupling, we can decompose the strong couplings 

as 

_ pQCD , loop _ pQCD(. , 

gj/i, = 9j/.^ +g.j/^ = gj/^ + 

g^, ^ g^?^^ + g'p^ ^ ± g^,) , (9) 

where g^^^^ and 5^?*^^ are couplings given by pQCD power counting, while gj°^ and g^?^ are given 
by an additional mechanism due to intermediate meson loop transitions; quantities 9,//,/, and q^i are the 
ratios of those two couplings for J/i/; and '0', respectively. Qualitatively, suppression of the g^i coupling 
implies that there exist large cancelations between g^'~^^ and 5^°^ while in J/?/; decays effects from 
may not be significant. 

Quantitative results supporting this require an explicit calculation of both g^fj^^ and g^'^^ , for which 
QCD models have been pursued in the literature. More or less, they respect the "12% rule" since they 
probe the charmonium wavefunctions at origin. The inclusion of intermediate meson loop transitions 
will introduce corrections to the couplings via the non- vanishing qj/^ and q^i in Eq. ([9]). It is worth 
noting that the couplings from the IML can be different for different decay channels. In particular, for 
pTT channel, it turns that |<z,7/^| < \qTp'\- 
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This relation again is not trivial at all. It further narrows down the meehanism that causes the 
deviations from the pQCD power counting, and also put a constraint on its behavior. As follows, instead 
of providing detailed calculations for the loops, we summarize the main features about the intermediate 
meson loop transitions and detailed numerical results will be reported later [3^ : 

I) Since both J/ip and ip' a-i'e below the open charm threshold, the intermediate meson loops will 
contribute to the real part of the couplings. This feature not only justifies the parametrization scheme 
in Step 2, but also makes the decomposition of the strong couplings in Eq. ^ physically meaningful. 

II) Since the '0' has a mass which is closer to the open DD threshold, its amplitude via the DD loop 
will be qualitatively larger than J/i/; due to near-threshold effects. 

III) Similar behavior due to intermediate DD{D*) and DD*{D) loops also shows up in a coherent 
study of J/ip and ip' jr/c and tp' — > 777^ [11] . 

IV) Light intermediate meson loops are strongly suppressed due to large off-shell effects. 

These features are consistent with a recent study of the "unquenched" effects arising from meson loops 
in Ref. ^33|, where it was shown that the intermediate meson loops still play an important role within 
charmonium states below the DD open threshold. 

V. SUMMARY 

In this proceeding, we carry out a systematic analysis of the problem of "p7r puzzle" , and clarify it on 
a more general ground. We show that the EM transitions play an important role in understanding the 
underlying mechanisms, which can be constrained by the isospin-violating channels. It thus allows us to 
separate out the EM amplitudes in those isospin-conserved channels. 

The nature of the VVP coupling as an antisymmetric tensor is also a key for disentangling the problem 
since whatever the mechanisms are for the transition, their contributions will simply be a correction to the 
coupling form factor. This allows a parametrization for the strong decay transitions in J/'ip{'ip') — » VP. 
The result shows that there exists an overall suppression on the strong decay amplitudes for ip' VP. 
Because of such a suppression, the strong decay amplitudes in some of those channels, such as pi: and 
K*K + c.c.^ become compatible with the EM transition amplitudes, with which the interferences produce 
significant deviations from the expectation of pQCD power counting rule. We then identify that the 
suppression on the strong decay amplitudes is originated from intermediate meson loop transitions, such 
as DD{D*), etc. In particular, the ip' is closer to the open DD threshold than J/4'- As a result, it will 
experience much larger threshold effects, and in this case a destructive interference. 

Such effects should be more general, hence may show up in other decay channels. We point out that 
a larger intermediate meson loop contribution originated from the same mechanism is also found in the 
study of J /ipiip') — > 7?7c and "0' m'c as an important mechanism interfering with the NRQCD leading 
amplitudes [331] . Nevertheless, the study of "unquenched" effects in the charmonium spectrum also gives 
rise to the importance of the intermediate meson loops [13]. 

In brief, we clarify that the "/ott puzzle" is not a single problem with the pi: channel. Instead, it 
is a rather general observation for all decay channels of J /ip{-4>') — > VP. In exclusive decays, multi- 
mechanisms can easily break down the pQCD power counting due to interference. In this sense, the 
so-called "p7r puzzle" should not be surprising. It turns to be more interesting to us that the systematics 
arising from a coherent study of all those VP decay channels can provide us with much deeper insights 
into the underlying dynamics. We expect that more precise measurement of those isospin-violating decay 
branching ratios at BESIII will help solve this long-standing problem [s^, . 
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